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ABSTRACT 

We study the constraints which the next generation of radio telescopes could place on 
the mass and number of neutrino species by studying the gravitational lensing of high 
rcdshift 21 cm emission in combination with wide-angle surveys of galaxy lensing. We 
use simple characterizations of reionization history and of proposed telescope designs 
to forecast the constraints and detectability threshold for neutrinos. It is found that the 
degeneracy between neutrino parameters and dark energy parameters is significantly 
reduced by incorporating 21 cm lensing. The combination of galaxy and 21 cm lensing 
could constrain the sum of the neutrino masses to within ~ 0.04 eV and the number 
of species to within ~ 0.1. This is an improvement of a factor of 2.6 in mass and 
1.3 in number over a galaxy lensing survey alone. This includes marginalizing over 
an 11 parameter cosmological model with a two parameter model for the dark energy 
equation of state. If the dark energy equation of state is held fixed at w = p/p = — 1 the 
constraints improve to ~ 0.03 eV and 0.04. These forecasted errors depend critically on 
the fraction of sky that can be surveyed in redshifted 21 cm emission (25% is assumed 
here) and the redshift of reionization (z = 7 is assumed here). It is also found that 
neutrinos with masses too small to be detected in the data could none the less cause 
a significant bias in the measured dark energy equation of state. 

Key words: large-scale structure of Universe - dark matter - gravitational lensing 
- intergalactic medium - low frequency radio astronomy 



1 INTRODUCTION 

One of the most tantalizing questions in experimental par- 
ticle physics - the nature of neutrino mass - is connected 
to one of the most tantalizing questions in observational 
cosmology - the nature of dark energy. Cosmological obser- 
vations presently put the m ost stringent upper limits on the 
absol ute mass of neutrinos (|Seliak et al . 2005; Spcrgc l~et al.l 
I2OO7I ) and might continue to do so for some time to come. 
However, a partial degeneracy between neutrino masses and 
the dark energy's equation of state limits how well either one 
can be determined from cosmological measurements based 
on the evolution of structure formation. 

Atmospheric and solar neutrino oscillation experiments 
strongly indicate that neutrinos have mass and that the sum 
of their masses is larger than the measured mass splittings 
5Z„ m v ~ \J Amj tm + yj Am'i un ~ 0.05 eV. Direct measure- 
ments from /3-decay experiments give an upper bound on 
the electron neutrino mass of 2.5 eV, th ere are no compet - 
ing measurements of the other flavors fsee lFogli et al.l (|2006l ) 
for a review of experimental results). Combinations of cos- 
mological data give a constraint of ^ v m v < 0.66 eV assum- 
ing a flat cosmology with a cosmological constant, but this 



limit is considerably looser if the density of da r k energy is 
allowed to evolve with ti me (|Seliak et al.l 120051 ; Irlannestadl 
l2005l ; ISpergel et alj|2007h . 

Massive neutrinos in the appropriate mass range are rel- 
ativistic when they decouple from the photons, electrons and 
baryons when the temperature of the universe is a few MeV. 
They have the effect of suppressing the power spectrum 
of density fluctuations on small scales (scales smaller than 
the horizon when the neutrinos become non-relativistic, 
k > fc nr ~ 0.018VmM m h Mpc -1 eV~ 1/2 where k is 
the Fourier wave number, m v is the neutrino's mass and 
h is the Hubble parameter in units of 100 km s~ Mpc -1 ) 
due to free-steaming. This suppression begins at the time 
of decoupling and continues until today because the ther- 
mal velocity of the neutrinos is still significant (v t hcim — 
150 (1 + z) m„ km s _1 eV _1 when non-relativistic). The free- 
streaming scale continues to decrease from its maximum of 



0.82 



h Mpc . Neutrinos 



(1 + zJ 2 "" V cV / 

will resist falling into dark matter halos with velocity dis- 
per sions ^ ^therm- For a r eview of neutrinos in cosmology 
see iLeseourgues fc PastorJ (|2006l ) . 

It has been proposed that neutrino masses could be 



2 R.B. Metcalf 



measured by the next generation of weak gravitational 
lensing surveys (EUCLIlflLSSTSPanSTA RR j| DE© 
(jHannestad et alj|2003 : iKitching et al.ll2008t ) and that the 
existence of massive neutrinos may limit the abil ity of these 
surveys to measure properti es of the dark energy (|Hannestadl 
120051 ; iKiakotou et al.|[2007h . Dark energy causes structure 
formation to evolve differently than it would otherwise at 
low redshift (z 1). A degeneracy arises between this and 
the late-time free-streaming of massive neutrinos since at 
scales significantly smaller than kfs the suppression is scale 
independent and the lensing surveys will be sensitive to a 
limited range in k. However, at z > 1 dark energy is ex- 
pected to have negligible effects on structure formation so 
if higher redshifts can be probed the degeneracy can be re- 
moved. The CMB (cosmic microwave background) provides 
information on the power spectrum at z ~ 1100, but Silk 
damping limits its sensitivity to scales below k nr . Gravita- 
tional lensing of pregalactic 21-cm radiation can provide a 
probe of structure formation on the scales needed, at the 
redshifts needed. 

The prospects for measuring gravitational lensing of the 
pregalactic 21 cm radiation have been s t udied by a num- 
ber of authors dZahn fc Zaldarriagal|2006l; 
2007bl;lHilbert et al.ll2007l;lLu fc Penll2007l; 



Metcalf fc White 



Metcalf fc White 



2007a! ) . The methods used in this pap er are described in 



more detail in 



aethods used m this pap 
iMetcalf fc White! (|2007ah 



2 FORMALISM & MODELING 
2.1 pregalactic 21 cm radiation 

A thoro ugh review of pregalac tic 21 cm radiation can be 
found in lFurlanetto et all |200^). For this paper we will men- 
tion only the barest essentials needed to specify our model 
for 21 cm emission. The fluctuations in the 21 cm bright- 
ness temperature depend on the spin temperature, T 3 , the 
ionization fraction, xh and the density of HI through 



ST b ~ 24(1 + 5 b )x H 



MB 



0.15 1 + z 



0.02 

1/2 

mK 



(1) 



n m h 2 10 

l|Fieldlll959l : iMadau et ai1ll997l ). As is commonly done, we 
will assume that the spin temperature is much greater than 
the CMB temperature. This leaves fluctuations in xh, and 
the baryon density St = (pt — 7>b)/7>b as the sources of 
brightness fluctuations. We will make the simplifying as- 
sumption that Xh = 1 until the universe is very rapidly 
and uniformly reionized at a redshift of «reion- We will 
take St to be distributed in the same way as dark matter 
accor ding to the CDM mode l. Nonlinear structure forma- 
tion (iPeacock fc Dodddll996l ) and linear redshift distortion 
l|Kaiser 1987 ) are included. Realistically, the reionization 
process will be inhomogeneous and may extend over a sig- 
nificant redshift range. This will increase C v (€) by perhaps 
a factor of 10 on scales larger than characteristic size of the 
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ionized bubbles (|Zaldarriaga et al1l2004l ) and make the dis- 
tribution non-Gaussian. 



2.2 gravitational lensing 

It is convenient to express the lensing formalism in terms of 
the convergence, k(8,z 3 ), at a position 9 on the sky which 
to an excellent approximation is related directly to the dis- 
tribution of matter through which the light passes 



k (ff, z^j = ^Ho 



" On f dz ^^g(z, Za )5 



with 



i 



/oo 
dz i) (z' , z 3 ) 



Zi+Sz 



dZ iJ E{zf g{z ' Zs) 
(3) 



D(z,0)D(z',z) 
D(z',0) ' 



(4) 



Ho is the Hubble parameter. The convergence can be 
thought of as a projected dimensionless surface density. The 
weighting function for the source distance distribution, 77(2), 
is normalized to unity. D(z' , z) is the angular size distance 
between the two redshifts and S(x, z) is the fractional den- 
sity fluctuation at redshift z and perpendicular position x. 
The function 



e{z) = ^£i m (i + zY + n A (i + zyf^), 



(5) 



where f2 m and f^A, are the present day densities of matter 
and dark energy measured in units of the critical density. 
It is assumed that the universe is flat - f2 m + = 1. The 
function describing the evolution of dark energy with red- 
shift can be written 



/(*) 



-1 



ln(l + z) 



[l + w(o)] din a (6) 

ln(l + z) 



where tu(a) is the equation of state parameter for the dark 
energy - the ratio of the of its pressure to its density - and 
a — (1 + z)~ is the scale parameter. 

For our purposes it is convenient to express equation Q 
as a matrix equation, 



K = G<5, 



(7) 



where the the components of K are the convergences run- 
ning over all position angles 9 and source redshifts, z 3 . The 
components of the vector S run over all position angles and 
foreground redshifts Zj. The matrix G is a function of most 
of the global cosmological parameters - Q m , Qa, w, etc - 
and is independent of position on the sky. This equation 
holds equally well if k and 5 are expressed spherical har- 
monic space or in the u-v plane where interferometer obser- 
vations are carried out. 

Within frequency bands of size 1 MHz we use an esti- 
mator for the Fourier transfo rm of the convergenc e , k(£, z s ), 
denoted k(£, z s ), based on the lZahn fc Zaldarriaga (2006) es- 
timator and then combine the bands as in IMetcalf fc White! 



Neutrino Masses, Dark Energy and the Gravitational Lensing of Pregalactic HI 3 



(2007a). The noise in k{£, v) within one frequency band is 
y, I d 2 £l [£-£'C v (£',k)+£-(£- 



c?(e',k)c?(\e!-e\,k) 



(8) 



where Cj (£, k) is the power spectrum of the actuai bright- 
ness temperature, while C u (£,k) = C v (£,k) + C„ (£,k) is 
the observed power spectrum which includes noise. £ is the 
angular Fourier mode number. Because the estimator is a 
sum over all the observed pairs of visibilities, it will (by the 
central limit theorem) be close to Gaussian distributed even 
though it is quadratic in the visibilities. 

The k(£, z)s can be grouped into a data vector 



D 



K 
K 



K 

GS 



(9) 
(10) 



where the components run over all the combinations of Zi 
and £ that are measured. To a good approximation modes 
from different bands and ^'s separated by more than the 
resolution of the telescope are statistically independent so 
the covariance matrix for D is diagonal 



Ny ^ 8ijN K (£i, Ui). 



(11) 



with \£i — £j\ larger than the resolution. The likelihood func- 
tion for the correlations in k(£, Zi) is given by 



ln/: = -iK t C- 1 K-i|C|, 



(12) 



where 



C = N + C K (13) 

l|Metcalf fc Whitell2007ah . C K here is the (cross-)power spec- 
trum of the convergence for two different source redshifts, 

[c] w = («(i; «*)«(!; *,■))• (14) 

This can be calculated using expression ([2]) and a model for 
the matter power spectrum. 

2.3 structure formation with massive neutrinos 

To calculate the lin ear matter power spect rum we use the 
analytic for mulae of lEisenstein fc Hul (|l999T ) with the mod- 
ification of iKiakotou et ah ( 2007 ) that improves accuracy 
w hen the neutrino masses are small. We use the method 
of iPeacock fc Doddsl (|l996l ) to transform this linear power 
spectrum into a power spectrum with nonlinear structure 
formation. This method has not been tested thoroughly 
against simulations of nonlinear structure formation with 
massive neutrinos, but we expect that the relation between 
the linear and nonlinear power spectrum will be only weakly 
affected by this. There is no baryon acoustic oscillations in 
the power spectrum. 

We assume that there are N u species of neutrinos with 
the same mass. The total density of neutrinos is fixed 
by the physics at the time of decoupling to be Q„h 2 = 

Wy/93.14 eV. If the masses are not degenerate the mea- 
sured N v will not be an integer. These are the parameters 
we will try to constrain. Roughly speaking N„ controls the 



average mass of the neutrino species and through this the 
free-streaming scale while fi„/i 2 affects the degree of sup- 
pression to the power spectrum. 

The measurements we consider are actually sensitive to 
any light particle species with significant cosmological densi- 
ties not just neutrinos. The interpretation of the parameters 
would be different for relic particles that were produced in a 
different way, for example axions, but the late-time physics 
would be the same. 



2.4 model observations 

2-4-1 21 cm observations 

We will concentrate on the planned SKA (the Square Kilo- 
meter Array jf] because it is the only planned telescope that 
will have a large enough size to be relevant to neutrino con- 
straints. It is only the core of the telescope that will be used 
for observing pregalactic 21 cm radiation. Plans for the SKA 
core have not been finalized, but it is expected to have a di- 
ameter of D tc i ~ 6 km (£ mstx ~ 10 4 ), an aperture covering 
fraction of / CO ver ~ 0.02 (the total collecting area of the 
telescopes divided by 7r(D te i/2) 2 ) and a frequency range ex- 
tending down to ~ 100 MHz which corresponds to z ~ 13. 
It is expected that the SKA will be able to map the 21 cm 
emission with a resolution of A8 ~ 1 arcmin. For reference, 
one arcminute (fwhm) corresponds to baselines of 5.8 km at 
z = 7 and 11 km at z = 15. We will take £ m i n = 10 to be 
the lowest mode to be measured. 

The noise in each visibility measurement will have a 
thermal component and a component resulting from imper- 
fect foreground subtraction. Here we model only the ther- 
mal component. If the telescopes in the array are uniformly 
distributed on the ground, the average integration time for 
each baseline will be the same and the power spectrum of 
the noise will be 



C 



2?r 



sys 



Avto \ /cover Dtcl J A^t / c 3 ovcr i' m ax(^) 2 ' 



(15) 



JZaldarriaga et all |2004| ; iMoralesI 120051 ; iMcQuinn et al.1 
|2006| ) where" T sys is the system temperature, Av is the 
bandwidth, t is the total observation time, and i? ma x(A) = 
2nD te \/X is the highest multipole that can be measured by 
the array, as set by the largest baselines. At the relevant 
frequencies, the overall system temperature is expected to 
be dominated by galactic synchrotron radiation. We will ap- 
proximate the brightness temperature of this foreground as 
T s k y = 180 K(i^/180 MHz) -2 ' 6 , as appropriate for region s 
well away from the Galactic Plane (|Furlanetto et aLll2006t ). 
This results in larger effective noise for higher redshift mea- 
surements of the 21 cm emission. We will consider an ob- 
servation time of 90 days which might be achievable within 
three seasons of observation and we will assume that the 
survey covers 25% of the sky. 

2-4-2 galaxy weak lensing survey 

For comparison and for combining with the 21 cm lensing 
we will consider a model galaxy lensing survey. The noise in 

5 www.skatelescope.org/ 
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power spectrum estimates from such a survey can be written 
as N K (£) = of /n e where n s is the angular number density of 
background galaxies and a e is the root-mean-square intrinsic 
ellipticity of those galaxies. This neglects all systematic er- 
rors as well as photometric redshift uncertainties. Following 
standard assumptions, we model the redshift distribution of 
usable galaxies as r/(z) oc z 2 e _( - z/ ' Zo ' , where z is set by 
the desired median redshift, and we adopt cr E = 0.25. The 
EUCLID satellite (the imaging part of which was previously 
known as DUNE0) proposes to survey 20,000 square degrees 
on the sky to a usable galaxy density of n g ~ 35 arcmin -2 
with a median redshift of z ~ 0.9. Several planned ground- 
based surveys - LSST, PanSTARRS - will cover comparable 
areas to EUCLID at a similar depth. In order to use tomo- 
graphic information, we divide the galaxies into 10 redshift 
bins each containing the same number of galaxies. 



2.4-3 CMB observations 

The Planck SurveyoiQ will do a full sky survey of the CMB 
radiation with higher resolution and more sensitivity to po- 
larization than is now available. We incorporate these future 
measurements into our forecasts by calculating the expected 
Fisher matrix. To do this we use the same te chnique as de- 
scribed in the appendix of lRassat et alj (|2008T l. This includes 
the contributions from both the scalar and tensor perturba- 
tions. 




Figure 1. The forecasted constraints on the combined mass and 
number of neutrinos with the fiducial model {J^ Tn v , Nv} = 
{0.66 cV, 3} marginalized over all the other cosmological param- 
eters. The dotted curves are for a EUCLID-like galaxy lensing 
survey, the dashed curves are for a 21 cm lensing survey, the solid 
curves are for the tomographic combination of the two surveys 
and the dot-dashed curve is for Planck by itself. 



3 RESULTS 

To assess how well the proposed observations could measure 
neutrino properties we adopt two statistical methods that 
are popular in the literature - Fisher matrix forecasts and 
the Bayesian evidence method. 



2.5 the cosmology 

We use all parameter cosmological model. The energy den- 
sities in dark energy and baryons are and Q.^ in units 
of the critical density. The density of matter is fixed to 
Q. m = 1 — Q.a to make the geometry flat. The primordial 
power spectrum is 



Po(k) = As 



H c- 



(16) 



The jfifj: parameter is included because if the primordial 
power spectrum is not a pure power-law it might partially 
mock the effect of early neutrino free-streaming on the power 
spectrum. The time dependent dark energy equation of state 
has two parameters w(a) = w + tu a (l — a). There are two 
neutrino parameter N v and X^„ m "- Including d "° and a 
non-constant w makes this a more general set of parameters 
with more potential degeneracies than has usually been used 
when predicting constraints on massive neutrinos. The op- 
tical depth to CMB last scattering surface is r. The fiducial 
model is set to {r, h, S!a/i 2 , flbh 2 , n s , -j^, w , w a , N u } = 
{0.09,0.7,0.343,0.0223,1,0,-1,0,3}. The normalization, 
A s , is set so that the fluctuations within a sphere of radius 
8 Mpc is erg = 0.75 in the fiducial model. Two values for 
S„m„ are used, 0.66 eV and 0.09 eV. All calculated param- 
eter constraints are marginalized over the other parameters. 



6 www.dunc-mission.net 

7 www.rssd.esa.int/indcx.php?project=Planck 



3.1 Fisher matrix forecast 

The maximum likelihood estimate for any parameter can be 
found by maximizing (|12p with respect to that parameter. 
The error in this estimator is often forecast using the Fisher 
matrix defined as 



d 2 \nC 

dpidpj 



(17) 



The expected error in the parameter p a , marginalized over 
all other parameters, is a 2 ~ (F -1 )^^. The unmarginalized 
error estimate (the error when all other parameters are held 
fixed) is (F aa )-\ 

Since £-modes separated by more than the resolution of 
the telescope will not be correlated, we can break the like- 
lihood function up into factors representin g each resolved 
region in £-space (|Metcalf fc White! l2007bl ). The result is 
that there are ~ (2£ + l)/ s k y independent measured modes 
for each value of £, where / s k y is the fraction of sky sur- 
veyed. The Fisher matrix can then be further simplified to 
the widely used form, 



= £ (2f+l)/ sky tr[C- 1 C, a C- 1 C, 



(18) 



The error estimates on the neutrino parameters are 
shown in figures [T] through [3] for combinations of different 
data and model assumptions. Figure[T]shows the constraints 
using the three kinds of data by themselves and the combina- 
tion of galaxy and 21 cm lensing. These lensing surveys are 
combined tomographically not just by adding their Fisher 
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3.0 

N„ 



Figure 2. The same as figure [T] except the Planck constraints 
are combined with the lensing constraints and the scale has 
been changed. There are two fiducial models at m v , N v } = 

{0.66 eV, 3} and m v , N v } = {0.09 eV, 3}. 



2.7 2.8 2.9 3.0 3.1 3.2 3.3 

N„ 



Figure 3. Same as figure [2] except the dark matter equation of 
state is held fixed at {w ,w a } = {—1,0}. 



S -1.00 




0.4 0.6 
£» m„ (eV) 



Figure 4. The constraints in the w - m v plane. Planck con- 
straints arc included. Note the improvement in the w a constraints 
especially when the neutrino masses are small. 



matrices as would be the case if they were independent mea- 
surements. Ten redshift bins with equal numbers of galaxies 
are used and ten redshift bins of fixed z-width are used for 
the 21 cm lensing. The CMB does a relatively good job of 
constraining the N v , but not as good a job of constraining 
the sum of the masses. 

In figure [2] the CMB constraints are combined with the 
lensing constraints. There is a drastic improvement because 
the parameter degeneracies are significantly reduced (Note 
the change in scale from figure [T}. The forecasted errors are 
0.04 eV and 0.1 for the combination of lensing surveys with 
m " = 0-66. This is an improvement of 2.6 in mass and 
1.3 in N v over either of the lensing surveys by themselves. 
In figure [3] the dark energy equation of state is kept fixed 
(w = —1). The improvement in the constraints illustrates 
the degeneracy between neutrinos and dark energy param- 
eters. The effect on the galaxy lensing constraints are par- 
ticularly strong. If one only considers cosmological constant 
models the constraints improve to 0.03 eV and 0.04 for the 
combined lensing case. 

Figure|4]shows the constraints in the w - ^ m " P iane - 
It can be seen more directly here how the degeneracy be- 
tween dark energy and massive neutrinos is reduced by in- 
cluding 21 cm lensing, especially when the neutrino masses 
are small. Note also that if it were assumed that neutrinos 
are massless, the maximum likelihood would not give the 
correct value for w . 

These calculations are in agreement with previous cal- 
culati ons of the constraints from future galaxy lens ing sur- 
veys (|Kitching et al.l 120081 ; lHannestad et al.l l2006t ) where 
more restricted cosmological models were used and thus 
stronger constraints were found. 



3.2 Bayesian evidence 

Another statistical question that could be asked is whether 
the data requires massive neutrinos. One way to answer this 
question is usi ng the average ratio of the Bayesian evidence 
lfjeffrevdll96ll ) for models with and without massive neutri- 
nos. This technique has be en used extensively i n cosmolog- 
ical p a rameter estima t ion (|Liddle et alj (|2006l ) ; ISaini et al.l 
|2004l ); lHeavens et all (1200717 for example) a nd galaxy lens- 



ing surveys specifically (jKitching et al 



12008ft . 



The Bayesian evidence is the probability of the data, 
D, given the model, M, marginalized over the parameters 
of that model, {#}, 



E(D\M)= / dO p(D\6,M)p(0\M). 



(19) 



The integral is over all of the parameter space. Bayes factor 
is the ratio of the evidences for two competing models 



B EE 



E(D\Mg) 
E(D\Mi) 



(20) 



Here model M is the simpler model with n parameters 
and model Mi is more complex and has more parameters, 
ni > n a . 

To make a forecast, Bayes factor must be averaged over 
expected data sets. It will be assumed that model M is the 
real case so that the averaging is according to this model. 
We are asking how well we can expect to role out model Mi . 
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In B 




0.2 0.3 

E m / eV 



Figure 5. The log of the average evidence ratio between a model 
with and without massive neutrinos. The axis are the deviation 
from the fiducial model {^ lv 'm v ,N l/ } = {0,3}. The three con- 
tour curves of each type are for 1, 2.5 and 5. According to Jeffreys 
(1961) |lnB| < 1 is "inconclusive", 1 < |lnB| < 2.5 is "substan- 
tial", 2.5 < | In _B| < 5 is "strong" and |lnB| > 5 is "decisive" 
evidence for the additional parameters. The expected constraints 
from Planck have been incorporated. The dotted curves are for a 
EUCLlD-like galaxy lensing survey, the dashed curves are for a 
21 cm lensing survey and the solid curves are for the tomographic 
combination of the two surveys. 



In B 




0.2 0.3 

E m / eV 



Figure 6. Same as figure \5\ except the dark matter equation of 
state is held fixed at {w ,w a } = {—1,0}. 



In calculating this we us e the approximation to (B) derived 
bv lHeavens et all (|2007f ) the Savage-Dickey ratio 



where 



66i 



F 



F 



(i); 



(21) 



(22) 



and the range of the indexes are i = l...n , a — n + l...ni 
and n,v = l...ni. We have assumed that there are no a 
priori limits on the range of the neutrino parameters. 

Figure [5] shows evidence ratio as a function of the devi- 
ation from the fiducial model, {X]„ m v , N„} = {0, 3}. Again 



we see that the constraints from the lensing of galaxies and 
21 cm put similar constraints on the neutrino parameters. 
By combining these data sets the constraints are improved 
by more than a factor of two in each dimension. A combined 
mass as small as 0.2 eV should be detected, 0.19 eV if N v 
is fixed to 3. Figure [6] shows how these constraints change 
if the dark energy equation of state is fixed to the cosmo- 
logical constant value, {w ,w a } — { — 1,0}. The constraints 
drastically improve which again demonstrates the degener- 
acy between dark energy and massive neutrinos. In this case, 
a combined mass of 0.09 eV would be detectable. 



4 DISCUSSION 

If it is assumed that neutrinos are massless, or have too 
small a mass to be of significance, the dark energy equation 
of state, w = p/p, could be underestimated from cosmologi- 
cal constraints. This would be the case even if the data does 
not allow for a detection of neutrino mass. It is also true that 
any independent constraint on the neutrino mass would im- 
prove future cosmological constraints on dark energy. The 
katrhnH beta-decay experiment expects to reach a level of 
0.2 eV for the electron neutrino mass and thus might have an 
impact on dark energy constraints. In section|3]it was shown 
that a neutrino mass as small as m v /N v — 0.03 eV could 
bias Wo high by ~ 1 a in future galaxy lensing surveys if it 
is not accounted for and if it is accounted for the error bars 
increase by a factor of 3.4 for w and 2 for w a . Lensing 
of high redshift 21 cm is a way to reduce this degeneracy 
by adding constraints on the growth of structure at higher 
redshifts where dark energy is presumed to contribute very 
little. Type la Supernovae surveys and surveys aimed at 
measuring the baryon acoustic oscillations at z ~ 1 will 
put constraints on dark energy that are based on the lumi- 
nosity distance redshift relation and so are independent of 
the growth of structure. It is important that methods based 
both on structure formation and on luminosity distance are 
fully realized to avoid systematic errors. The two probes are 
also necessary to test alternative theories of gravity on large 
scales which might also provide an explanation for the ap- 
parent acceleration of the cosmological expansion. 
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